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I. INTRODUCTION 



Gyromagnetic ratios of particles in hydrogenlike bound states have become in the last 
ten-fifteen years an active field of experimental and theoretical research. The gyromagnetic 
ratio of a bound electron is proportional to the ratio of the spin flip and cyclotron frequencies 
of a hydrogenlike ion and to the electron-ion mass ratio. The experimental uncertainties of 
the ratio of the spin flip and cyclotron frequencies of the hydrogenlike carbon 12 C 5+ and 
oxygen 16 7+ were reduced to 5 — 7 parts in 10 10 , see [H, @] arid review in The theoretical 
expression for the bound state ^-factor was also greatly improved recently (see, e.g. 0, 0] 
and references in [3|), and the theoretical uncertainty was reduced to 1.5 — 5.5 parts in 
10 11 . As a result measurements of the bound electron g-factor became the best source for 
precise values of the electron mass in atomic units [3( • This bright picture is marred by the 
discrepancy on the magnitude of the leading relativistic corrections to bound state ^-factors 
existing in the literature. Leading relativistic (both nonrecoil and recoil) corrections to g- 
factors of loosely bound spin one half particles were calculated a long time ago @-@] ■ These 
corrections in the case of loosely bound particles with arbitrary spins were calculated in j^] . 
It turned out that the leading relativistic and recoil corrections to bound state ^-factors are 
universal, they do not depend on the spin of the constituents but only on their free ^-factors. 
Among other results, this feature allowed the authors of [9] to sum all nonrecoil and recoil 
corrections of order a n (Za) 2 to bound state ^-factors in hydrogenlike ions. Universality of 
leadin g bi nding corrections to bound state ^-factors was challenged in [l(J The results 
of [l(J 11] contained the contributions explicitly depending on spin, both for nonrecoil and 
recoil corrections. These terms shift the theoretical value of the bound state ^-factors of the 
hydrogenlike carbon 12 C 5+ and oxygen 16 7+ by about 2 — 3 parts in 10 11 . The discrepancy 
will become even more phenomenologically relevant if proposed improvement [l2| of the 
experimental accuracy bytwo orders of magnitude is achieved. Later the universal results 
of 0] were confirmed in 13. . 



Below we derive an effective nonrelativistic QED (NRQED) Hamiltonian for charged par- 
ticles with arbitrary spins and calculate the leading relativistic and recoil corrections to the 
bound state ^-factors in loosely bound two-particle systems. We show that these corrections 
are universal for all spins, and explain the physical reasons behind this universality. 



II. NRQED LAGRANGIAN 

A loosely bound two-particle system is effectively nonrelativistic, with characteristic ve- 
locities of constituents of order Za. We are looking for the leading nonrecoil and recoil 
corrections to bound state ^-factors of order (Za) 2 . NRQED is a natural framework for 
calculation of these corrections. The NRQED Lagrangian sufficient for calculation of these 
corrections should include terms in nonrelativistic expansion up to and including terms of 
order v 2 . The well known (see, e.g., [uj ]) NRQED Lagrangian for the spin one half case 
is constructed from the covariant derivatives do + ieA and D = V — ieA = i(p — eA), 
electric and magnetic fields E and B, and the spin operator S. Notice that in a loosely 
bound two-particle system the scalar potentialA is of order v 2 , {eA ) ~ (Za) 2 . For higher 
spin particles we should, besides the spin, include also higher irreducible intrinsic multipole 
moments as the building blocks of the NRQED Lagrangian. Technically these multipole mo- 
ments are polynomials in the components of the spin operator that for higher spins do not 
reduce to numerical tensors and operators linear in spin. The general NRQED Lagrangian 
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has the form (compare with 15] for spin one half case 
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where Qy = SiSj + SjSi — (2/3) S 2 Sij is proportional to the electric quadrupole moment op- 
erator (Qij = for spin one half), and is a 2s + 1-component spinor field for a particle with 
spin s. We included in the Lagrangian in Eq. ([T]) operators of dimensions not higher than 
four, but omitted some terms (like the terms with derivatives of magnetic field3) that are 
irrelevant for calculation of the leading recoil corrections. Let us mention that gauge invari- 
ant bilinears in E and B are of too high order to generate leading relativistic contributions 
of order {Za) 2 to bound state ^-factors. 

The coefficients in the NRQED Lagrangian for spin one half charged particles are de- 
termined from comparison of the scattering amplitudes in nonrelativistic theory and in 
relativistic QED. We would like to follow the same path for the case of arbitrary spin, but 
renormalizable QED for charged particles with high spin does not exist. However, the rules 
for calculation of all one-photon interactions of charged particles with arbitrary spin were 

Trj ] . This construction uses only Lorentz invariance and 



constructed some time ago in [16| 
local current conservation, and is valid for charged particles of arbitrary spin. The interac 
tion vertex that includes all higher multipole moments in the approach of [if 
generalization of the ordinary spin one half vertex 



17 is a direct 
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2m 2m 
where q = P2 — Pi, is the generalization of ordinary spin one half a uu , 
spin four-vector, r — (q ■ S) 2 , and F e (0, 0) = 1, F m (0,0) = g/2. The wave functions are 
spinors with dotted and undotted indices that are symmetrized among themselves (for more 



Sfj, is a covariant 



details see 16|, |17J, |l9j)- The form of the vertex in Eq. ([2]) is uniquely fixed by the require- 
ments of Lorentz invariance, C, P and charge conservation. Charged particles with higher 
spins automatically carry higher multipole moments that arise as coefficients in expansion of 
the form factors F e and F m over r [l6, 17j]. These intrinsic electric and magnetic multipole 
moments are treated phenomenologically, and we do not try to calculate them. The phe- 
nomenological approach to multipole moments is an advantage for our purposes because we 
would like to describe how ^-factors of not necessarily electromagnetic origin (for example 
the g-factor of a spin one deuteron) change in a loosely bound electrodynamic system. 
We find coefficients in Eq. ((Tj) by comparing one-photon scattering amplitudes in NRQED 



and in the relativistic formalism of [16|, LL7J with vertex Eq. (J2]). Although some terms in 



Eq. ([T]) are bilinear in electromagnetic fields A and E they still can be restored from one- 
photon terms due to gauge invariance. Therefore the one-photon relativistic vertex in Eq. 



1 See [l3| for the explicit form of the Hamiltonian that includes all operators with dimensions not higher 
than four. 
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is sufficient for calculation of all the coefficients in Eq. ([I]). We calculated scattering ampli- 
tudes off an external electromagnetic field using the nonrelativistic Lagrangian in Eq. (TjQ) 
and using the relativistic one-photon vertex in Eq. (j2J) at q 2 = and r = 0. In the rel- 
ativistic calculation we used noncovariantly normalized particle spinors in the generalized 
standard representation, which is necessary for consistency with the respective nonrelativis- 
tic results. Diagrammatically this choice of spinors and representation corresponds to the 
Foldy-Wouthuysen transformation (for more details, see, e.g., [l9j). After nonrelativistic 
expansion we compared results of the relativistic calculation with the nonrelativistic ones 
and obtained values of all constants in the Lagrangian in Eq. (TjQ) 



g £2 
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where £ 2 = 4S, X = 1/(25 - 1) for integer spin and S 2 = 45* + 1, A = 1/(25) for half 
integer spin. Dependence on the magnitude of charged particle spin arose in the coefficients 
before the Darwin term and the induced electric quadrupole interaction. The ^-factor in 
Eq. ([3]) is the gyromagnetic ratio defined by the magnetic form factor F m (0, 0) = g/2 in 
Eq. (TJ]). In the spin one half case g reduces to a sum of the QED perturbation series 
if the charged particle is subject only to electromagnetic interactions. The coefficients in 
Eq. ([3]) in the spin one half case coincide with the respective coefficients in [15|, if the 
phenomenological g-factor is substituted in the expressions in 1151 instead of the perturbative 



g = 2(1 + a/27r). As an independent test of the effective Lagrangian in Eq. ([I]) we considered 
the charged W^ ± -boson sector of the Glashow-Weinberg-Salam Electroweak Theory amended 
by the anomalous magnetic moment term. We derived the effective NRQED Lagrangian for 
the W ± bosons. This Lagrangian coincides with the Lagrangian in Eq. ([T]) for spin one 
charged particles. 

The coefficients in Eq. fl3]) are calculated ignoring all loop diagrams and q 2 and r de- 
pendence of the form factors in Eq. fl2]). Both the loop diagrams in relativistic QED and 
multipole expansion of the form factors would generate further corrections to the coefficients 
in Eq. 03]). In ordinary renormalizable spin one half QED (as well as in the renormalizable 
QED of spin one vector bosons) all diagrams, besides those that give contributions 
only to the free particle g-factors, generate corrections to the coefficients that are addition- 
ally suppressed by powers of Za. We expect the same effect in any reasonable theory for 
higher spin particles. It is also obvious that accounting for q 2 and r dependent terms in 
the form factors in Eq. ([2]) generates terms suppressed by additional powers of Za. No- 
tice that we ignored the internal electric quadrupole moment implicit in the form factor 
F e (q 2 ,r) in Eq. (j5J), but still a term with an induced electric quadrupole moment arose in 
Eq. (JTJ) with a coefficient dependent on the magnitude of spin. This immediately means 
that electric quadrupole interaction with external electric field depends on the magnitude of 
spin. Something like this mechanism could in principle make the interaction of the magnetic 
dipole moment with an external magnetic field spin-dependent, and lead to dependence of 
bound state g-factor on the magnitude of spin. But this happens neither with the coefficient 
before the term S B, nor with the coefficients before the other terms in Eq. ([1]) that give 
contributions to leading binding correction to the bound state ^-factors. 



4 



Formally two-photon relativistic Compton effect diagrams are needed to obtain coeffi- 
cients before all terms in the NRQED Lagrangian bilinear in electromagnetic fields. How- 
ever, as mentioned above, all gauge noninvariant terms in Eq. (JT|) bilinear A and E can be 
restored from one-photon diagrams with the help of gauge invariance. Any gauge invariant 
terms connected with the two-photon diagrams are of too high order in Za to contribute to 
the leading relativistic corrections of order (Za) 2 . This allows us to avoid consideration of 
the relativistic two-photon diagrams. 

The Lagrangian in Eq. ([1]) with the coefficients from Eq. ([3]) is sufficient for calculation 
of the leading relativistic corrections to the bound g-factor in the nonrecoil case. 



III. EFFECTIVE TWO-PARTICLE HAMILTONIAN. CENTER OF MASS MO- 
TION 

Next we construct the effective nonrelativistic quantum mechanical Hamiltonian for a 
loosely bound electrodynamic system of two particles needed to calculate both nonrecoil 
and recoil corrections of order (Za) 2 to the bound state g-factors. The interaction between 
two charged particles with accuracy up to (Za) 2 is described by the one photon exchange 
which generates Coulomb and Breit interactions. We calculated the two-particle scattering 
amplitude for particles with arbitrary spins and magnetic moments using the vertices from 
Eq. (J2J) and the photon propagator in the Coulom b g auge. After nonrelativistic expansion 



we obtained the interaction potential (see also [13l . Il8l 
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where T"i( 2 ), Pi(2), ^1(2), ^1(2) > 9i(2), an d are the coordinate, momentum, spin, mass, 

gyromagnetic ratio, and induced quadrupole moment of the first (second) particle, and 
r = 7"i — r 2 is the relative coordinate. 

This interaction is a natural generalization of the spin one half one-photon potential (see, 
e.g., [HI). The only difference is that like in the Lagrangian in Eq. ([T|) the coefficients 
in Eq. (j4j) before the Darwin terms depend on the magnitude of particles' spins, and new 
terms with induced electric quadrupole moments arise (we again ignore intrinsic electric 
quadrupole moments and other intrinsic multipole moments not relevant for calculation of 
the leading corrections to bound state g factors). The interaction in Eq. (j3J is calculated 
in the absence of a small uniform external magnetic field that is present in the g-factor 
problem. This drawback is easily repaired by the minimal substitution pi — > pi — e^Ai, 
A l = Bx n/2. 

Combining the nonrecoil Lagrangian in Eq. ([1]) and the one-photon potential (after mini- 
mal substitution) in Eq. (j4j) we obtain a total effective nonrelativistic two-particle quantum 
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mechanical Hamiltonian for electromagnetically interacting particles with arbitrary spins 
(we preserve below only the terms relevant for calculation of the ^-factor contributions) 



where 



H = H 1 + H 2 + H int , (5) 
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and if 2 is obtained from Hi by the substitution 1 — > 2. 

The nonrelativistic effective two-particle Hamiltonian in Eq. (JSJ) describes all (nonrecoil 
and recoil) leading relativistic corrections to bound state ^-factors of each of the constituents. 
The coefficients before all terms in the Hamiltonian in Eq. (jSJ) do not depend on the magni- 
tude of spin, and the corrections to the bound state ^-factors are connected with the terms 
in the Hamiltonian that are linear in external magnetic field. Already at this stage we see 
that these corrections are universal and do not depend on the magnitude of spin. 

Actual calculation of the recoil corrections requires separation of the effects connected 
with the motion of the bound system as a whole from the internal effects. This task is not 
quite trivial because the center of mass variables do not separate in the presence of external 
field. For the current case of a small external magnetic field a solution was suggested in 
0, il]. The main idea is to impose the condition that the center mass of a loosely bound 
system moves in a small external field exactly in the same way as a respective elementary 
particle with the same mass and charge. Neither canonical nor kinetic momentum of a 
charged particle are conserved in external uniform magnetic field. Instead, in the symmetric 
gauge A(r) = B x r/2, the pseudomomentum p + qA(r) is conserved [6|, 

[H,p + qA(r)] = 0, (8) 

where p and r are the canonical momentum and coordinate of the charged particle, q is its 
charged, and the Hamiltonian has the standard form H = (p — qA{r)) 2 / (2m). Classically 
conservation of pseudomomentum means that the center of the Larmor orbit remains at 
rest. A naive transition to the standard center of mass coordinates r = T\ — r 2 , R = 
H\ri + /i 2 r 2 (/ij = mi/{mi + m 2 )) does not secure conservation of total pseudomomentum 
P + ( e i + e 2)A(-R), where P = pi + p 2 . To satisfy the transparent physical requirement 
of total pseudomomentum conservation (conservation of the position of the Larmor orbit 
center) transition to the center of mass coordinates should be accompanied by the unitary 



All charges in this paper carry sign, so, for example, for an electron e 
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transformation of the Hamiltonian H — > U 1 HU, where U = e i(eiM2 e wi) A ( R )- r , After this 
transformation the internal Hamiltonian acquires the form 

H = + ^ + <l + H?X + H r , (9) 

where m r = mim 2 / (mi + m 2 ) is the reduced mass, H^ n (H^ n ) describes interaction of the 
first (second) spin with the external field, and H r includes all other terms in the Hamiltonian. 
Explicitly the Hamiltonian for the first particle spin interaction with the external magnetic 
field is 
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,(10) 



The Hamiltonian i^ sp - n for the second constituent spin has a similar form. The leading 
binding correction to the (/-factor is completely described by the Hamiltonian in Eq. (110p . 
We calculate its matrix element with the help of the first order perturbation theory between 
the Schrodinger-Coulomb wave functions that are eigenf unctions of the unperturbed internal 
Hamiltonian. After simple calculations we obtain the bound state (/-factors with account of 
the leading relativistic corrections of order {Za) 2 for s-states with the principal quantum 
number n 
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These results resolve the discrepancy mentioned in the Introduction in favor of the results 
in (9[ (see also jUJ). The remarkable property of the expressions in Eq. (fill and Eq. ( IT2l 
is that they are universal for particles of any spin; they depend only on the (/-factors of free 
charged particle, not on the magnitude of their spins. Technically this happened because no 
terms in the effective two-particle NRQED Hamiltonian in Eq. (JSJ) relevant for calculation 
of the leading relativistic corrections contain spin- dependent coefficients Aj, E 2 . 
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IV. UNIVERSALITY AND THE BARGMANN-MICHEL-TELEGDI EQUATION 



Universality of the (Za) 2 corrections to the bound state g-factors in Eq. ( ITT]) and Eq. (Fl2|) 
requires physical explanation. We already noticed above that simple analysis of the dimen- 
sions and spin structure of all terms in the NRQED Lagrangian in Eq. ([T|) leads to the 
conclusion that terms with derivatives of electric fields do not generate contributions to the 
leading relativistic corrections to the bound state g-factors. Omission of the field deriva- 
tives is the basic assumption for validity of the Bargmann-Michel-Telegdi (BMT) equation 



19|,|20j. Hence, one can use the BMT equation for motion of spin in external electromagnetic 
field for derivation of the leading corrections to the bound state g-factors as was suggested in 
[9| . Then the leading relativistic corrections are universal if the BMT equation is universal 
for all spins. To explain universality of the coefficients in the BMT equation let us recall 
the main steps in its derivation. The basic idea is to use Lorentz invariance to generalize 
the nonrelativistic precession equation 

S = 9 ^SxB (13) 
to a relativistically invariant equation for motion of spin in external field. Following 19[ 2 



we introduce a four-pseudovector a M (a 2 = — S 2 ) to describe spin. In the particle rest frame 
a M = (0, S). In the relativistic generalization of Eq. ( TT3j) the derivative on the LHS should 
be over proper time dr = ^dt (7 = 1/Vl — v 2 ). On the RHS we are looking for a vector 
that is linear and homogenous in the gauge invariant external EM field F^ v and in the spin 
pseudovector a^. Then the most general relativistically invariant equation for spin motion 
is 

— = aF" v a v + fiv?F vX u v a x , (14) 
dr 

where w M is the four-velocity, = (j,jv), u^u^ = 1, and a and /3 are unknown constants. 
The values of these constants are uniquely restored using the nonrelativistic precession 
equation Eq. ffT3]) and the classical relativistic equation of motion for a charged particle 
in external field. We obtain 

,-±-a~ to -2)£. (15) 

Returning to ordinary time t and noncovariant spin vector S we can write the BMT equation 
in the form 



We see that all coefficients in the BMT equation are universal for any spin, and this univer- 
sality follows from the universality of nonrelativistic spin precession and Lorentz invariance. 

Next we sketch the derivation of the leading relativistic corrections to bound state g- 
factors based on the BMT equation [3], in order to demonstrate how the universality of the 
BMT equation leads to universality of the leading corrections to the bound state ^-factors. 
At the first step we represent Eq. ffTol) as the Heisenberg equation of motion for the spin 
operator 



S 



with the Hamiltonian 



(17) 
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To calculate the leading relativistic corrections to the bound state ^-factor we expand 
the Hamiltonian in Eq. (|T8|) up to quadratic terms in u 2 , 7 « 1 + v 2 /2. We obtain the 
nonrelativistic Hamiltonian 

fl ,-A(f 9 >-f) s . £ - (9 - 2) ((p-^- B )( s '(p-^ 



2m [ \ m? J 2m? 

S .^x( P -e A )] | 
m J 

where we made the substitution = (p — e A) / m required by gauge invariance. The Hamil- 
tonian in Eq. ( TT9|) coincides with the one in Eq. (Q , and explains why the latter Hamiltonian 
has universal coefficients. 

Besides the Hamiltonian in Eq. ([H]) calculation of the recoil corrections to bound state g- 
factors also requires knowledge of the interaction Hamiltonian in Eq. ([7]). Above we derived 
this interaction Hamiltonian from the relativistic one-photon exchange, but it is easy to see 
that it includes only the ordinary nonrelativistic spin-orbit and spin-other orbit interactions. 
It is well known that these interactions are universal for any spin and we could write the 
Hamiltonian in Eq. ([7]) using only nonrelativistic quantum mechanical expressions. This 
finally explains universality of the corrections of order {Za) 2 to the bound state g-factors. 
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